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Abstract 

Graham and Sloane proposed in 1980 a conjecture stating that ev- 
ery tree has a harmonious labelling, a graph labelling closely related to 
additive base. Very limited results on this conjecture are known. In 
this paper, we proposed a computational approach to this conjecture by 
checking trees with limited size. With a hybrid algorithm, we are able to 
show that every tree with at most 31 nodes is graceful, extending the best 
previous result in this direction. 

1 Introduction 

In 1980, Graham and Sloane first introduced the notion of harmonious labelling 
on general graphs in [5]. For a general graph G — (V, E) with v nodes and e > v 
edges, a harmonious labelling for G is an injection / : V — >• Z e such that each 
edge x,y G E receives a distinct sum f(x) + f(y) in Z e . For trees, as we have 
e = v — 1, we demand / to be a surjection instead. If G allows a harmonious 
labelling, we say that G is harmonious. 

Graham and Sloane proposed the following conjecture. 

Conjecture 1. Every tree is harmonious, that is, for every tree T = (V,E) 
with n = Card(V), there is a surjection f : V —¥ % n -i such that g : E — » 
Z„_i, {x, y} i-> f(x) + f(y) is a bijection. 

According to Gallian's dynamic survey on graph labelling [4] , there are few 
results on harmonious trees. Graham and Sloane proved that caterpillars are 
harmonious. Aldred and McKay proved in jT] that every tree with no more than 
26 nodes is harmonious by computing a harmonious labelling for each such tree. 

In this paper, we use a similar computational approach to prove the following 
result. 

Theorem 1. Every tree with at most 31 nodes is harmonious. 

We extend the result in [1] by proposing and implementing a new algorithm 
for harmonious labelling. This algorithm is very similar with the one in a recent 
work of the author [5] on graceful labelling that also improved previous result. 
We will discuss details of this algorithm in the following section, which will be 
followed by implementation details, observations and discussions. 

In the following, we will use T = (V, E) to denote a general tree, with n = 
Card(V) denoting its number of nodes. We call any surjection f : V —> 7L n -\ 
a labelling of T, and for any v £ V we call f(v) the label of v. We will call the 
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accompany g : E — > Z n _i, {x,y} i-> /(a;) + /(y) the induced labelling of /, and 
17(e) the edge la&e/ of e € £. Operations on labels are in Z„_i unless otherwise 
stated. 

2 Algorithm 

To find harmonious labellings, we use a hybrid algorithm consisting of three 
sub-algorithms, each searches in different ways. The first two sub-algorithms, 
probabilistic backtracking and tabu search, resemble their counter parts in the 
previous work of the author [3]. The third, two-stage constraint solving, is new. 

All three sub- algorithms have the property that they may fail to return a 
labelling, but when they succeed, the returned labelling is harmonious. This is 
in fact favourable, since we do not know whether every tree is harmonious, and 
if there is one in our range, our algorithm should not be blocked. 

We use the method in [6] to generate free trees with a certain size, thus trees 
are represented by its level sequence. Trees are rooted in one of their center, 
which is midpoint of any longest path. Level of each node is its distance to the 
chosen root, and level sequence is a canonical sequence of levels of nodes. We 
can identify nodes with terms in the sequence. 

2.1 Probabilistic Backtracking 

By backtracking, we mean that we try to assign labels to nodes successively in 
some order, and when a node has no valid label, we backtrack to previous nodes 
for alternative choices. 

We choose to label nodes in the order of level sequence, which has a good 
property that every node is linked to some node before in the sequence. When 
we label nodes in this order, each assignment of label fixes a certain edge label. 
We label the root arbitrarily. 

When we try to label a node, we may have various choices, but some of them 
may violate the condition for harmonious labelling. We say that a label is valid 
for the current node to be labelled if it maintains the property of the partial 
labelling on nodes other than the root being an injection, and also the partial 
induced labelling. It is easy to see that, if we always assign valid labels in each 
step, we will obtain a harmonious labelling. 

We set up empirically a limit on the number of backtracking to avoid long 
runtime due to combinatorial explosion. This limit fails the algorithm when ex- 
ceeded. To avoid being trapped in local minimal, random restart, which restarts 
the algorithm after a certain iterations, and stochastic pertubation, which ran- 
domly swaps pairs from time to time, are also used in actual implementation. 

Here is the pseudocode of this algorithm in Algorithm [T] 

2.2 Tabu Search 

Finding harmonious labelling can be modelled as a combinatorial optimization 
problem. For a labelling /, which is already a surjection, we define the following 
evaluation function. 

Eval(f) = n - 1 - Card({g(e)\e 6 E}) 
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Algorithm 1 Backtracking Search 
Assign random label to the root 
while iteration limit is not exceeded do 
if Current node has valid label then 

Pick one at random to assign and move to the next node 

else 

Backtrack to previous node 
end if 

if all nodes are labelled then 

return the obtained labelling 

end if 
end while 
return failure 



We have Eval(f) = exactly when the induced labelling g is an bijection, 
that is to say / is harmonious. We want to find a labelling / that gives as 
low value of Eval(f) as possible, which turns out to be a typical combinatorial 
optimization problem. This is also the approach in [1]. 

We use tabu search to solve this problem. We start from a randomly chosen 
labelling /. The replicated label does not matter, as for any harmonious la- 
belling /, / + c with c G Z n _i is still a harmonious labelling. Then we attempt 
to lower the value of Eval(f) by swapping labels between nodes. If a swap 
between u, v can lower Eval(f), we keep the change and put {it, v} into a tabu 
list, which forbids this swap to be done again in several turns. This is where 
"tabu" comes from. 

We also stop computing when the number of iterations exceeds a certain 
limit, with the same reason to avoid long runtime. 

Here is a pseudocode in Algorithm 



Algorithm 2 Tabu Search 

Choose an arbitrary labelling / which is already a surjection 
while iteration limit is not exceeded do 

Randomly pick several pairs of nodes not forbidden 

if some pair u,v such that Eval(f) decreases by swapping f(u),f(v) 

then 

Swap f(u) and f(v) 

Forbid u, v from being chosen in a few next iterations with tabu list 
end if 

if / is a harmonious labelling then 
return / 

end if 
end while 
return failure 
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2.3 Two-stage Constraint Solving 

Finding harmonious labelling can also be modelled as a constraint satisfaction 
problem (CSP) on variables f(v),v G V, with the following constraints: 

f:V^{0,l,...,n-l} 

Alldiff (f(v),vcV) 

Alldiff(f(u) + f{v) mod (n - 1), {u.v} e E). 

Alldiff stands for "all different", which means its terms should be distinct. 
Clearly this model expresses exactly the condition of harmonious labelling, mod- 
ulo a constant. Our previous backtracking algorithm reduces to a routine ran- 
domized backtracking algorithm to solve CSP in this model. More information 
about CSP solving algorithms can be found in a book by Apt [2J. 

However, in this form, it is not easy to implement some usual solving tech- 
niques, such as forward checking or look ahead. We try to reduce the problem 
to a nicer form by the following two-stage solving. 

Let L be the set of leaves of T, and we note P(u) the parent of some leaf u. 
In the first stage, we find a partial labelling / over V — L that does not violate 
conditions of harmonious labelling. This can be solved with the backtracking 
algorithm, since we have similar situation. With this partial labelling, we can 
reduce the model to the following form. 

/':L^{0,l,...,n-l}-7m(/) 

Alldiff (f'(v),v 6 L) 
Vv G L, f'(v) + f(P(v)) mod (n -l)£G, f'(v) $ Im(f) 

where G = {g(e)\e € E, e n L = 0} the set of existing induced labels. 

This model decouples largely dependences between all /'(«)• They are now 
linked only by an Alldiff constraint. In this model, we can effectively perform 
the forward checking optimization, which shrinks domains (possible values of a 
variable) by removing impossible values after each fixation of a variable. This 
is not effective in the first model, because variables share the same domain and 
the valid labels of a node depends on the value of its parent. 

However, the new model is also limited by the chosen partial labelling. Even 
if there is a harmonious labelling, we cannot find it unless we guessed its partial 
labelling over V — L correctly. To solve this problem, we repeat this algorithm 
for several times. At each run, if we solved the reduced model, we obtain 
a harmonious labelling by extending the partial one. If we could not find a 
harmonious labelling in a certain number of runs, we simply return a failure. 

2.4 Hybrid Algorithm 

The principle idea of this hybrid algorithm is that we filter trees through each 
sub- algorithms. Only trees failing the current sub-algorithm are sent to the next. 
The order is two-stage constraint solving, backtracking, then tabu search. This 
order is determined empirically with decreasing speed and increasing possibility 
of finding a harmonious labelling. 

If a tree fails all three sub-algorithms, it may be a good candidate of counter- 
example of the conjecture of Graham and Sloane. 
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3 Implementation and Result 



The previous hybrid algorithm is implemented in C++, including tree generation. 
All parameters are fixed empirically. There are various changes of algorithms 
for optimization, but they are all minor. Using this implementation, we are able 
to check computationally whether some finite families of trees are harmonious. 

On an Intel Core 2 Duo T7200, we checked that every tree with at most 
31 nodes is harmonious. We found a harmonious labelling for every tree with 
our hybrid algorithm. Verification of trees with exactly 31 nodes took roughly 
2800 hours of CPU time. To obtain the result in [T] , only 8 hours of CPU time 
is needed. 

Empirically, runtime for a single tree grows exponentially with its size, but 
only slightly and tolerable comparing to the growth of number of trees of the 
same size. With more computational power, we can extend further our checking 
range. 
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